SM and MSSM Higgs Boson Production: Spectra at large transverse Momentum 
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Strategies for Higgs boson searches require the knowledge of the total production cross section and the trans- 
verse momentum spectrum. The large transverse momentum spectrum of the Higgs boson produced in gluon 
fusion can be quite different in the Standard Model and the Minimal Supersymmetric Standard Model. In this 
paper we present a comparison of the Higgs transverse momentum spectrum obtained using the Pythia event 
generator and the HiGLU program as well as the program Hqt, which includes NLO corrections and a soft 
gluon resummation for the region of small transverse momenta. While the shapes of the spectra are similar for 
the Standard Model, significant differences are observed in the spectra of Minimal Supersymmetric Standard 
Model benchmark scenarios with large tan l3. 

PACS numbers: 14.80.Bn, 14.80.Cp, 12.60.Jv 



I. INTRODUCTION 

The search for Higgs bosons belongs to the most impor- 
tant endeavors at the Large Hadron Collider (LHC) in order 
to establish experimentally the Higgs mechanism for elec- 
troweak symmetry breaking. In the Standard Model (SM) one 
isospin Higgs doublet is introduced, which leads to the exis- 
tence of one physical Higgs particle after electroweak sym- 
metry breaking, while the other three degrees of freedom are 
absorbed by the W and Z bosons fr|. In the SM Higgs sec- 
tor the only unknown parameter is the Higgs mass. Based on 
triviality and unitarity arguments its value is required to range 
below ~ 800 GcV |2, 3]. The Higgs couplings to fermions 
and electroweak gauge bosons grow with the corresponding 
masses (V = W, Z) 
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where Gp = 1.16637 x IQ-^ GeV~^ is the Fermi con- 
stant. Therefore, the Higgs couplings to the W and Z bosons 
as well as to third-generation fermions are phenomenologi- 
cally relevant, while the couplings to the first two genera- 
tions are less important. The direct search at the LEP2 ex- 
periments excluded Higgs masses below 114.4 GeV |4]. If 
the SM is embedded in a Grand Unified Theory (GUT), the 
quadratically divergent radiative corrections to the Higgs self- 
energy tend to push the Higgs mass towards the GUT scale 
Mgvt ~ 10^^ GeV. In order to establish a Higgs mass of 
the order of the electroweak scale an unnatural fine-tuning of 
the counter terms is required. This hierarchy problem remains 
unsolved within the SM. 

The most attractive solution to the hierarchy problem is 
the introduction of supersymmetry (SUSY), a novel symme- 
try between fermionic and bosonic degrees of freedom |3l- 
Due to the additional contributions of the SUSY partners of 
each SM particle the quadratic divergences in the Higgs self- 
energy are canceled. The hierarchy problem is solved ja], 
if the SUSY particle masses are maximally of the order of 
a few TeV. In the minimal supersymmetric extension of the 
SM (MSSM) two isospin Higgs doublets have to be intro- 
duced in order to preserve SUSY |7] and to render the model 



TABLE I: Higgs couplings in the MSSM to fermions and gauge 
bosons (V — W, Z) relative to the SM couplings. 
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free of anomalies. After electroweak symmetry breaking five 
Higgs bosons are left as physical particles: two CT'-even neu- 
tral (scalar) particles h, H, one CV-odd neutral (pseudoscalar) 
particle A and two charged bosons H^. 

At leading order the Higgs sector is determined by two in- 
dependent input parameters, which are usually chosen as the 
pseudoscalar Higgs mass Ma and tan/? = f2/i'i, the ra- 
tio of the two vacuum expectation values. The light scalar 
Higgs boson h has to be lighter than the Z boson at lead- 
ing order. This upper bound, however, is significantly en- 
hanced to a value of ^ 140 GeV due to radiative correc- 
tions, which are dominated by top- and stop-loop contribu- 
tions LL^. Moreover, all Higgs couplings are affected by 
the same type of coiTections. The couplings of the Higgs 
bosons to fermions and gauge bosons are modified by coeffi- 
cients, which depend on the angles a and j3, where a denotes 
the mixing angle of the two CT'-even Higgs fields. The cou- 
plings, normalized to the SM Higgs coupling, are listed in Ta- 
ble|I] An important property of these couplings is that for large 
values of tan/3 the down(up)-type Yukawa couplings are 
strongly enhanced (suppressed). The direct Higgs searches 
at the LEP2 experiments have excluded neutral Higgs masses 
MhM^^'^-^GeY and Ma< 91.0 GeV as well as charged 
Higgs masses Mh± < 78.9 GeV 1 10]. 

At the LHC the dominant neutral SM (MSSM) Higgs pro- 
duction mechanisms (for small and moderate values of tan /3 
in the MSSM) are the gluon fusion processes I 111 

gg^h {h, H, A) 




FIG. 1 : Leading order contribution to the SM process gg — > h. 



which are mediated by top and bottom loops (see Fig. [Q and 
for the neutral CT'-even MSSM Higgs bosons h, H in addition 
by stop and sbottom loops with the latter contributing signif- 
icantly if the squark masses are below ^ 400 GeV |12]. At 
large values of tan /3 Higgs radiation off bottom quarks be- 
comes competitive within the MSSM ilj. The NLO QCD 
connections to the top and bottom loops enhance the cross 
sections by 50-140% for the SM Higgs boson 1 14, 15] and 
the MSSM Higgs particles for small values of tan/3 fllfT^l . 
while for large values of tan (3, where the bottom loops pro- 
vide the dominant contributions due to the strongly enhanced 
bottom Yukawa couplings (see Fig.|2ji, the corrections amount 
to 10-60%. It should be emphasized that the NLO QCD cor- 
rections are of more moderate size, if the bottom loops be- 
come dominant, while they are large in the regions of top- 
loop dominance. The NNLO coiTections are known in the 
heavy top mass limit, which is a valid approximation only for 
regions, where the top loops are dominant, i.e., for the SM 
and the MSSM for small values of tan/?. In this heavy top 
mass limit the top loop reduces to an effective ggh coupling 
(see Fig.|3jl.' This limit is expected to provide a good approx- 
imation to the exact total cross section if Mh ^ 2mt. In the 
SM the maximal deviation from the fully massive NLO result 
is less than ^ 5% in this mass range, while in the MSSM it 
increases to ~ 30% for tan /? ~ 5 fnlfl^ . if the full mass de- 
pendence of the LO cross section is taken into account, while 
the X-factor is derived in the heavy top mass limit. Compared 
to the LO result, the cross section is enhanced by a factor of 
1.7 to 2.3 at NNLO fT?]. Moreover, the SUSY-QCD coiTec- 
tions have been calculated within the MSSM in the limit of 
heavy SUSY particle and top masses. They turn out to be large 
for the squark loops 1 12], while the genuine SUSY-QCD cor- 
rections, mediated by virtual gluino and stop exchange, are 
0(5%) and thus small ]20]. In the gluon fusion processes 
the Higgs bosons are produced with vanishing transverse mo- 
menta at leading order. 

Higgs boson production at finite transverse momenta re- 
quires the additional radiation of a gluon or quark in the gluon 
fusion process. The leading order contributions to the differ- 
ential gluon fusion cross section da/dpx stem from diagrams 
as those in Fig. |3^ Analytical results can be found in 12 1|] . 



' Note that due to the top Yukawa couphng the top contribution does not 
decouple for large top masses, but approaches a mass-independent value. 

^ In this letter we do not consider Higgs radiation off bottom quarks 
gg, qq — > bbH, which is of comparable size as the gluon fusion process 
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FIG. 2: Dependence of the K factors for the gluon-fusion cross sec- 
tions on the value of tan /3. The corresponding K factors obtained 
by omitting the bottom loops are: Kh = 1.71, Kh ~ 1.76 (Mr = 
150 GeV), Ka = 1.78 {Ma = 150 GeV), Kh = 1.91 {Mh = 
500 GeV) and Ka = 1.87 {Ma = 500 Ge^ independent of 
tan/3. CTEQ6L1 (CTEQ6M) parton densities (ij are used for the 
LO (NLO) cross sections with the corresponding Higgs mass as the 
renormahzation and factorization scale. 
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FIG. 3: Effective ggh coupling in the heavy top mass limit. 



NLO contributions to these expressions are again only known 
in the heavy top mass limit f22], which serves as an approx- 
imation for Mfi,PT ^ "T-t- Since the bottom mass is small, 
this approximation is not expected to work in MSSM regions, 
where the bottom loop contributions are significant. At low 
transverse momentum (pr <C Mh) multiple soft gluon emis- 
sion spoils the validity of fixed order calculations. To ob- 
tain reliable results in this region, contributions of all orders 
have to be taken into account by resumming the large loga- 
rithmic terms ln"'{m\/p'^). These computations have been 
performed at leading logarithmic (LL), next-to-leading log- 
arithmic (NLL) 1 23] and next-to-next-to-leading logarithmic 
(NNLL) level f24]. The NNLL calculation was performed 
in the heavy top mass limit. To obtain reliable predictions 
for the differential cross section over the whole pr range, the 



for large values of tan /3 within the MSSM. 
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FIG. 4: LO Feynman diagrams contributing to the transverse momentum spectrum of the SM Higgs boson, mediated by Higgs couplings to 
gluons. 



resummed results are matched to fixed order calculations at 
large transverse momentum. This is for example done in the 
Hqt 125II program (for more details on the programs see be- 
low), results of which are shown in Fig.|5] The heavy top mass 
limit is only valid as long as pT ^ mt- For larger values this 
approximation tends to overestimate the cross section f2l|]. 
This can be seen in Fig.|6j where the SM spectrum was com- 
puted with HiGLU l26ll for the exact top mass and an infinite 
top mass. 



have the same shape in the SM and the MSSM. A crosscheck 
of the results of Pythia at large transverse momenta is per- 
formed with the program HiGLU [26]. In the SM the spectrum 
obtained by HiGLU is softer than that of Pythia, because 
Pythia uses the heavy top mass approximation. The assump- 
tion in Pythia of similar transverse momentum shapes in the 
SM and the MSSM is not correct in general, especially in re- 
gions of bottom-loop dominance 1 29|. In some MSSM bench- 
mark scenarios large discrepancies between the two programs 
are found. 
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FIG. 5: Comparison of the NLO differential cross section with the 
NNLL result, matched to the NLO result for SM Higgs bosons with 
mass Mh = 115 GeV in the heavy top mass limit. The renormaliza- 
tion and factorization scales are chosen as the transverse mass. 



The development of strategies for Higgs boson searches at 
the LHC requires reliable estimates of the total and the differ- 
ential cross section. The latter can potentially be used to dis- 
criminate signal from background, especially at large trans- 
verse momenta, where many background processes are sup- 
pressed [27]. In this letter we compare the transverse mo- 
mentum distributions of the SM and MSSM Higgs bosons 
obtained with the programs PYTHIA |28], HiGLU |26|] and 
Hqt II25I] . In Pythia the large transverse momentum spectra 



II. NUMERICAL RESULTS AND DISCUSSION 

If not mentioned otherwise, all results are given for a Higgs 
boson mass of 115 GeV and a value of tan/3 = 30. For the 
numerical studies the following programs have been used: 



Pythia v6.227 [28] 

The Pythia process gg -^ h uses the matrix element at 
order a^, i.e., the Higgs boson is produced at rest. The 
Higgs boson receives its transverse momentum only by 
initial state radiation, which is added to the hard inter- 
action. The initial state radiation is added in such a 
way that the large pT spectrum is matched to that of 
the process gg — > hg. All calculations implemented in 
Pythia are performed in the heavy top mass limit. We 
have used the default scale choices. 



- HiGLU |2e 

HiGLU is a program to compute the total Higgs pro- 
duction cross section via gluon fusion at NLO in the 
SM as well as in the MSSM. It can also be used to cal- 
culate the differential cross section at the same order 
in Og (i.e. at LO), although this option is not docu- 
mented. To the best of our knowledge, HiGLU is the 
only freely available program which allows to com- 
pute the differential spectrum in the MSSM with the 
full heavy quark mass dependence. The renormaliza- 
tion and factori zatio n scales are chosen as the transverse 
mass Mt = \/mI 
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FIG. 6: Comparison of the differential spectrum in the SM generated 
by HiGLU and Pythia for Mh = 115 GeV. The renormalization 
and factorization scales in HiGLU are chosen as the transverse mass. 



FIG. 7: Pythia spectra for the SM and the four benchmark scenarios 
scaled to the SM total cross section. The light scalar Higgs spectra 
are shown for A/^ = 115 GeV. 



Hqt|25|1 

Hqt can be used to compute the px spectrum of the 
Higgs boson produced in gluon fusion at LO and NLO. 
The spectra are calculated in the heavy top mass limit, 
but they can be normalized to the exact total cross sec- 
tions. Hqt only performs calculations in the SM. It in- 
cludes a soft gluon resummation, thus providing a reli- 
able and finite prediction in the limit of small transverse 
momentum in contrast to the purely perturbative result 
implemented in HiGLU, which diverges for pr — > 0. 
The renormalization and factorization scales are chosen 
as the transverse mass Mt — \f^^ 



Pt- 



If not mentioned otherwise, CTEQ5L parton densities were 
used in PYTHIA, while in HiGLU and HQT CTEQ6M parton 
densities were used L17il . 



A. Standard Model 

As shown in Fig. |5] good agreement between the spectra of 
Pythia and Higlu is obtained in the SM, if Higlu is used 
in the heavy top mass limit. The shapes of the two curves 
coincide reasonably well at large pT- At low pT the HiGLU 
result diverges, because HiGLU does not perform any resum- 
mation. While good agreement is achieved in the heavy top 
mass limit, significant deviations for large values of pT are 
observed when taking into account finite quark mass effects 
in Higlu: the cross section at large pt is overestimated by 
the heavy top mass approximation 1 2 1 ] as can be inferred also 
from Fig.|6l Thus any analysis requiring Higgs bosons at large 
transverse momenta has to be adjusted accordingly. 



B. Minimal Supersymmetric Standard Model 

In the MSSM the comparison between Pythia and Higlu 
was carried out for the four different benchmark scenarios 



proposed in |30], i.e., the rin}^'^^, the no-mixing, the gluopho- 
bic Higgs and the small tteff scenario. 

In Pythia the shape of the spectra is the same for the 
SM and all MSSM scenarios, only the normalization changes. 
This is true even though by switching on supersymmetry in 
Pythia, the couplings of the quarks to the Higgs boson 
change compared to the SM. As in the SM, the Higgs particle 
is produced at rest and then initial state showers are added to 
the hard scattering. This shower is using a splitting kernel that 
is convoluted with the LO gg -^ h matrix element, thereby 
ensuring that the shower reproduces the gg -^ gh matrix ele- 
ment at large px values by a proper matching procedure. This 
convolution of the splitting kernel with the hard interaction is 
performed in the heavy top mass limit in the SM and MSSM. 
Therefore the spectrum shapes coincide in both cases. This 
can be clearly seen in Fig. where all spectra produced by 
Pythia, scaled to the corresponding total SM cross sections, 
are shown. 

In contrast to this, HiGLU shows differently hard spectra 
for the benchmark scenarios, as can be inferred from Fig. |S] 
All spectra turn out to be softer than the SM spectrum. Tak- 
ing the Pythia spectrum as the base of a physics analysis, 
the Higgs signal can be overestimated by more than an order 
of magnitude for px > 100 GeV. The softness of the spectra 
can be traced back to the fact, that the bottom quarks yield 
the main contribution to the differential cross section for large 
tan j3. Tableimiists the couplings for the four benchmark sce- 
narios. In the m™"^ scenario, the bottom loop contribution 
is enhanced by a factor of 1.65 x lO"* compared to the top- 
loop contribution alone. A lighter quark in the loop generates 
a softer spectrum as can be read off Fig. |5] The same effect is 
also visible in the spectra of the other scenarios. It is particu- 
larly large in the small a^s scenario. 

The effect of two other parameters, the Higgs mass and 
tan/3, on the large pT spectrum is shown in Figs. 1101 and 
[TTI If the mass of the light Higgs boson approaches its up- 
per limit (which corresponds to Ma -^ oo), the spectrum 
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FIG. 8: Comparison of the differential spectra generated by HiGLU and PYTHIA in the four benchmark scenarios for the light scalar Higgs 
boson with A'lh = 115 GeV. 



TABLE II: LO cross sections and Higgs Yukawa couplings to up- 
and down-type quarks in the four benchmark scenarios computed 
with Higlu for tan /3 = 30 using CTEQ6L1 parton densities for the 
light scalar MSSM Higgs boson with Mh = 115 GeV. 
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The transverse momentum distributions changes by ^^ ±40%. 
This variation can be taken as a lower bound of the theoret- 
ical uncertainty at LO in analogy to the total inclusive gluon 
fusion cross section. 



C. Best estimate for the large pr spectrum 

The aim of this subsection is to get the best possible predic- 
tion for the differential cross section at large transverse mo- 
mentum. It is clear from the preceding discussion that the 
following conditions have to be met: 



becomes SM-like. The same happens if tan/3 is lowered as 
presented in Fig.fH] since in both limits. Ma -^ oo and small 
tan /3, the Higgs Yukawa couplings become similar to the SM 
couplings. Fig. [21 shows the variation of the spectra if the 
renormalization and the factorization scales are v aried by a 
factor two around the transverse mass Mt ~ V-^'^h + Pr- 



• The full quark mass dependence has to be taken into 
account 

• For large tan/3 the bottom quark loops must not be 
omitted 

• The calculation has to be performed at the highest pos- 
sible order in a,. 



Higlu m =~ 
Higlum = 174.3 GeV ^ 
Higlu m = 100 GeV 
HIglu m = rrijj 




400 500 

P^ „. [GeV] 

■^T, Higgs 



> 


10 


a> 




O 


1 




n 




a. 








H 


10' 


Q. 




■a 




B 


10-^ 


■a 






10-^ 




10-^ 



10 = 



10' 



, 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 = 

Higlu SM E 


i\ 


Higlu mj^'", tan p = 3 = 


- \k^ 


Higlum"''", tan |3= 10: 


^ \^ >- 


Higlu m™', tan p = 20 1 


^\ '~^~^ 


Higlu m"'^^ tan p - 30 - 


= ^s<. 


->;::^„..__^ Higlu m|^^^ tan p = 50 = 


^so- 


"'-■C^- 


= > 


v^"-. "■^'^■"■-^ = 












^>^^""-. '"'■■' — ■■■-^ 


E 


■^>0*--.. ■■••-^:;"*^ 






" 


^^>-\- -. ^ 


= 


, 1 , , , , 1 , , , !^^>^ 



100 



200 



300 



400 



500 



P^ „. [GeV] 

■^T, Hrggs "■ 



FIG. 9: Differential cross section in the SM for different top masses 
andMh = 115 GeV. 



FIG. 11: Differential light scalar Higgs cross section for different 
values of tan /3 and Alh = 115 GeV. All spectra are reseated by the 
ratio of the total LO SM cross section and the corresponding total 
LO cross section of the scenario. 
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FIG. 10: Differential light scalar MSSM Higgs cross sections for 
different Higgs masses in the in™"^ scenario. 



By using HiGLU, the two first conditions are automatically 
met. The third one is not completely fulfilled, because HiGLU 
performs the calculation only at third order in the strong cou- 
pling constant. Presently there are no calculations at order a* 
fulfilling the first condition, since at this order the differential 
cross section is only known in the heavy top mass limit. In the 
meantime, the following improvement to the HiGLU result is 
the only possibility. From the comparison of the differential 
cross section in the heavy top mass limit at LO and NLO a px 
dependent A'-factor can be extracted, which then can be ap- 
plied to the exact LO result. This procedure has been carried 
out for the example of the m™"^ scenario. The pT dependent 
iiT-factor was computed with Hqt and then applied to the LO 
Higlu spectrum. The result is shown in Fig.^] The figure 
shows the Hqt spectra at LO and NLO, from which the K- 
factor was extracted, the HiGLU result at LO as well as the 
scaled HiGLU spectrum. For comparison the Pythia result 
is also given. However, it should be noted that this approxi- 
mation will be valid only for small values of tan (3, where the 
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FIG. 12: Light scalar transverse momentum spectra in the m™"^ sce- 
nario for different renormalization (fi) and factorization (Q) scales 
md Mh = 115 GeV. 



top loops are dominant, and not too large transverse momenta, 
while its validity for large values of tan /?, where the bottom 
loops become dominant, is not clear, before a fully massive 
NLO calculation is available for large transverse momenta. 



III. CONCLUSION 

The large transverse momentum spectrum of the Higgs bo- 
son produced in gluon fusion was investigated with different 
programs. While the predictions of Pythia and HiGLU agree 
reasonably well within the SM, significant differences were 
found in the MSSM at large values of tan/3. Compared to 
Higlu, Pythia overestimates the differential cross section 
in some scenarios by more than one order of magnitude for 
Pt > 100 GeV. The reason for this is, that Pythia is always 



Scaled Higlu 

Higlu MSSMm = 174.3G6V 

HQT SM LO m, = «. 
HQTSMNLOm, =~ 
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500 



working in the heavy top mass limit, where the quark loop 
is replaced by an effective coupling of the two gluons to the 
Higgs boson. In this way, the shape of the Pythia spectrum 
is not sensitive to changes of the Higgs Yukawa couplings, 
which modify the relative weight of the top and bottom quark 
loop contributions. 
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